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1 Introduction 

The initial transverse energies and multiplicities in central (zero impact parameter) 
ultrarelativistic heavy ion collisions have been computed in |l|] from perturbative QCD 
supplemented by the crucial assumption of saturation of produced semihard gluons: a 
simple saturation criterion defines a saturation scale psat ~ 1 (2) GeV at RHIC (LHC). 
Doing the computation at transverse momenta pt > Psat gives an estimate of the 
effect from all transverse momentum scales, both above and below Psat- "Initial" here 
then means proper times of the order of 1/psat ~ 0.2 fm (RHIC) and ~ 0.1 fm (LHC). 
Assuming thermalisation at formation and further entropy conserving expansion, these 
initial gluon numbers can be converted to final hadron multiplicities. The predicted 
multiplicities agree well with the first results from RHIC 0. 

The results of ^ are formulated in the form of scaling rules, quantity ~ CA°'{^)^, in 
which the constants C, a, h are determined for central AA collisions. A further variable 
one has experimental control over is the centrality or impact parameter dependence of 
multiplicities and transverse energies. In fact, in |0] the centrality dependence of the 
multiplicity has been proposed as a means of distinguishing between various models 
for particle production in ultrarelativistic heavy ion collisions. In particular, in one 
observed that the ratio of charged particle multiplicity to the number of participants, 
Nch/Npa^T-t, grows when b decreases while it decreases in the HIJING Monte Carlo model 
1^. The purpose of this note is a systematic study of the behaviour of charged particle 
multiplicities with impact parameter in the pQCD + saturation model. 

2 Impact parameter dependent saturation criterion 

Since the present work is a simple extension of |l| to non-central b 7^ collisions, we 
refer there for a detailed exposition of dynamical ideas. 

The tools of describing non-central AA collisions are [Q the nuclear density n^(r'), 
normalised to A, the standard nuclear thickness function Tyi(6) = / dznAir), normalised 
as J (PbTA{b) = A, and the nuclear overlap function T^Aib), J d?bTAAip) = A"^, which 
is a 2-dimensional convolution in b-space of two T4(&)'s. We shall systematically use 
the physically relevant Woods-Saxon nuclear density distribution with Ra = 1.12 A^^^~ 
0.86 y4~^/^ fm, for which these integrals have to be computed numerically. 

Experimentally, however, one does not have a direct control over the impact param- 
eter. Instead, using forward calorimeters one can measure the number of participants 
A^part in the collision; for central collisions A^part^2y4. This can be approximately con- 
verted to impact parameter by the formula 
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(6) = J dhTAih - s) [1 - exp(-ai,rB(s))] + J dhTB{s) [1 - exp(-ai„T^(b - s))] 
= 2j d^sTA{h-s)[l-expi-a,^TA{s))] {A = B), (1) 



where cxin ~ 40 mb is the inelastic pp cross section in the RHIC energy range. An 
example is shown in Fig. ^ The Poissonian exp(— cTinT^) in ([^) can as well be replaced 
by a binomial-like (1 — aT /A)^] the numerical effect is negligible. We have taken 
(Jin = 35, 39 and 42 mb for y/s = 56, 130 and 200 GeV, respectively, based on 0. 

We have now the tools to generalise the b = saturation condition for initially 
produced partons, 

iV^A(Po,b = 0) = 2T^^(0)(TpQCD(Po)= E Mpo)=pIR\, (2) 

k=g,q,q 

to arbitrary b. Perturbative QCD enters via an LO computation of the inclusive 
production cross section <JpQCT){po) of (mini)jets with pr > po and \y\ < 0.5, with NLO 
contributions taken into account via an overall i^-factor K = 2 according to ; all the 
formulas are spelled out in detail in [^. Eq. (^ expresses the fact that at saturation 
Naa{Po = Psnt) quanta each with transverse area vr/p^g^^ fill the whole nuclear transverse 
area 7tR\. Numerical, group theory factors or powers of g could be included 0, but 
these are anyway 0{1) unless one discusses a parametric weak coupling limit g ^ 0. 
All parton flavours are included, though at po = psat gluons clearly dominate even at 
lowest energies. 

In terms of an average transverse area density in central collisions, dN^A/d'^s ~ 
NAA{Po^O)/{7rR\), Eq. (||) becomes dNAA/d'^s ■ tt/pI = 1. For arbitrary b and s, the 
average saturation criterion thus generalizes to a local one, 

dNAAiPoMs) ^ 2:r^^b - s)TAis)a,^cBiPo) = pl/n. (3) 

We shall choose the coordinate system at b 7^ so that b = (b, 0) and so that the 
centers of the nuclei are at (6/2,0) and (—6/2,0). A solution of Eq.(||) gives us the 
local saturation momentum Psat{b, s) at the position s of an A+A collision at impact 
parameter b. Fig. |l| shows Psa.t{b = 0,s), which only depends on s; Fig. |^ shows 

Psat(6 = Ra, S = {Sx, 0)) and Psat(6 = Ra, S = (0, Sy)). 

The predicted initial multiplicity at fixed b then simply is dNAA/d'^s at the solution 
of (^ integrated over s. In terms of the right-hand side: 

NAA{b) = I d'spUb, s)/7r. (4) 

Here an important issue enters: how does one deal with large values of s or small values 
oi Psa.t{b,s) in (^? Physically, how does one treat very peripheral collisions for which 
Psat becomes small, nonperturbative? One may note that even at SPS energies, where 
Psat ~ 0.7 GeV, the data was reproduced by this model up to a 20% error. Clearly the 
pQCD+saturation model can be reliable only if the large s, small psat region makes a 
negligible contribution to the integral (|^). This, in fact, is the case: in central collisions, 
for example, the contribution of the range from psa.t > 0.7 GeV to psat > 0.5 GeV to 
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Figure 1: The saturation nionientuni Psat(& 
b = 0, ^ = 197 and ^ = 56, 130, 200 GeV. 



0,s) computed as a solution of Eq.(^) for 



the integration (|) is only an Ri 4% increase in Naa- In Figs. |T| and ^ the decrease in 
p1^^ for psat < 0.5 GeV is very rapid and the surface integral (Q) is negligibly affected. 
Values of N^Aip) as computed from (H) including only p^ax > 0.5 GeV are plotted in 
Fig. ^. The total initial particle production becomes thus described by the saturation 
model in an effective way in the sense that by pushing the local saturation criterion to 
its limits (to psat -^ 0) no additional, non-saturated, components need to be included. 
For very peripheral collisions, or for pp collisions, the model cannot be applied. 

The values of h in Figs. |l| to ^ can be converted to the experimentally more directly 
accessible number A^part using Eq. ([1|) and Fig. |^. Concretely, Figs. [^ to |^ extend to 
h = I.SRa, which corresponds to A^part ~ 30; similarly b = I.GRa would correspond to 
A'part ~ 50. For still larger b one expects saturation effects to disappear and this region 
is thus beyond the scope of the current study. 

We shall here not discuss the transverse energy Et production in detail, since for it 
the relation between the initially produced and finally observed Et per unit rapidity is 
likely to be more complicated than that for multiplicity, due to pdV work and transverse 
flow effects in the course of expansion. However, a reasonable estimate for the initial 
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Figure 2: The saturation momentum Psat{b = Ra, s) along the symmetry axes of the overlap 
area. Solid curves are for s = {sx,0) and the dotted ones for s = (0, Sy), computed from Eq. 
(^ for a Woods-Saxon nuclear density distribution for A = 197 and ^/s = 56, 130, 200 GeV. 

production is the same Et per quantum as at b = 0: ET{h)/N{h) ^ 1.35 (psat(^)), 
where (psat(&)) is computed with a weight TA(s)Tyi(b — s) for Psat(&, s). Note also that 
it is here that these pQCD+saturation results deviate significantly from the classical 
field computations [0, |Tl|: their Et/N is larger by a factor 3. 



3 Comparison with experiment 

The above were predictions for the multiplicity integrated along a curve Tj ~ l/psat(&, s) 
of initial time. Comparing the multiplicities with experiment requires assumptions on 
the expansion and decoupling of the system. These questions are discussed extensively 
in the literature and one is looking forward to experimental tests. 

We shall here use the simplest assumption of early thermalisation and entropy con- 
serving longitudinally boost invariant expansion. At the central slice z = the lo- 
cal particle density is initially then n{Ti{b, s),s) ^ dN/cPs/^TiArj) = 5(rj,s)/3.6, 
where rj is the space-time rapidity and S is the entropy density. The initial to- 
tal entropy per rapidity unit. Si, can thus be computed from Si{b) = 3.6A^aa(&)- 
This predicts that the relation of the number of initially produced particles (glu- 
ons) to the total number of particles in the final state (pions per unit rapidity) is 
Nf/NAA = {Nf/Sf)/{NAA/Si) * {Sf/Si) ^ 3.6/4, inserting the number/entropy ra- 
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Figure 3: The initial (r w 1/psat) number of quanta NAA{b) produced in ^^^Au+^^^Au 
collisions at ^/s =56, 130 and 200 GeV as a function of b/RA (the upper set of three curves), 
as computed from Eq.(Q) with psat ^ 0.5 GeV. The three lower curves show A'part(^) computed 
from (H) for cjin = 35,39,42 mb (bottom to top). The average values of A^part(&) given by 
PHOBOS H are indicated by the two (almost overlapping) arrows on the vertical axis, the 
values of A^aa(&) at the corresponding values of h/RA are marked by the asterisks. 



tios for massless gluons and massive pions. To obtain the experimentally measured 
dNch/dr], ?7=pseudorapidity, we use the estimates A^ch ~ |^/ and dN/dr] ^ 0.9dN/dy. 
Fig. ^ then shows the quantity {dNch{b) / drj) / {Npi^j.t/2) , i.e., the number of charged 
particles produced in unit pseudorapidity per number of participant pairs, as a func- 
tion of the number of participants. The parameters {A and ^/s) are those for RHIC 
experiments. 

For central collisions (approximately; as shown in Fig.|^, the number of participants 
in is somewhat less than 2A) the agreement with experiment is good. The rest of 
the curve is a prediction. 

The prediction in Fig.^has the striking feature of being essentially constant in A'^part- 
This is in agreement with the estimate in 0, shown by the dotted lines. The small 
difference is due to the fact that in |^ one did not actually perform the saturation 
computation at arbitrary b but instead used b = results at Aes = A^part(b)/2. Due 
to the computed ^ A-scaling of TVaa at b = 0, A^aa ~ A"-^^^, this immediately results 
in A^AA/^part ~ l/^part' ^ ^low iucrcasc as A^part decreases. At very small A^part this 
would seem to lead to a striking effect, but then one enters the very peripheral region 
in which saturation does not work. 
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Figure 4: Solid curves: the final number of charged particles per unit pseudorapidity and 
per participant pair produced in Au+Au collisions at y^ = 56, 130, 200 GeV vs the number 
of participants. Dotted curves: the same for central collisions with effective nuclei ^gfj = 
A^part/2. The PHOBOS data for 56 and 130 GeV are also shown, with separate error bars 
for statistical and systematic errors [Qj. 

4 Conclusions 



In this note we have extended the pQCD + saturation model for initial production of 
partons in A + A collisions, studied in quantitative detail in [|l|, from central (b = 0) 
to non-central collisions, up to perhaps b = 1.6Ra. This extension is based on a 
simple geometric argument. Subject to the usual uncertainties in relating initial to 
final observed numbers, predictions for particle multiplicities and transverse energies 
can now be made and tested in terms of three variables. A, ^/s and b (or number of 
participants). 

It is impossible to give a systematic estimate of the theoretical error in the pre- 
dictions. Corrections to pQCD can, in principle, be computed, but the use of the 
saturation scale adds a nonperturbative, even phenomenological, element. What one 
can say is that the model is unambiguously determined, fits the data at SPS within 
20% and predicted well the first measurements at RHIC, so maybe it is an acceptable 
way of at least extrapolating to so far unobserved parameter values. It thus offers a 
strongly constrained framework for an analysis of experimental data in a 3-parameter 
{A,y/s, b) space and one should be able to make a distinction between different models 
along the lines suggested in |Q. 
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